The beryllium atom and beryllium positive ion in strong magnetic fields 
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The ground and a few excited states of the beryllium atom 
in external uniform magnetic fields are calculated by means 
of our 2D mesh Hartree-Fock method for field strengths rang- 
ing from zero up to 2.35 • 10 T. With changing field strength 
the ground state of the Be atom undergoes three transitions 
involving four different electronic configurations which belong 
to three groups with different spin projections Sz ~ 0, —1, —2. 
For weak fields the ground state configuration arises from the 
ls^2s^, Sz — configuration. With increasing field strength 
the ground state evolves into the two Sz — —1 configurations 
ls^2s2p_i and ls^2p_i3(i_2, followed by the fully spin po- 
larised Sz ~ —2 configuration ls2p_i3d_24/_3. The latter 
configuration forms the ground state of the beryllium atom 
in the high field regime 7 > 4.567. The analogous calculations 
for the Be^ ion provide the sequence of the three following 
ground state configurations: ls^2s and ls^2p_i (Sz = —1/2) 
and ls2p_i3d-2 {Sz = -3/2). 

I. INTRODUCTION 

The behaviour and properties of atoms in strong mag- 
netic fields is a subject of increasing interest. This is 
motivated by the astrophysical discovery of strong fields 
on white dwarfs and neutron stars ||l|-^. On the other 
hand the competition of the diamagnetic and Coulomb 
interaction causes a rich variety of complex properties 
which are, of course, also of interest on their own. 

For a long time the investigations in the literature fo- 
cused on the hydrogen atom (for a list of references see, 
for example, |Q-0|). As a result of the corresponding in- 
vestigations the absorption features of certain magnetic 
white dwarfs could be understood in detail and a mod- 
elling of their atmospheres was possible (see ref. Q for a 
review up to 1994 and M for more recent references). De- 
tailed spectroscopic calculations were carried out recently 



for the helium atom in strong magnetic fields pd . These 
calculations allow to identify spectra of other, namely 
helium-rich objects, including the prominent white dwarf 
GD229 ||ll|. Recently a number of new magnetic white 
dwarfs have been found whose spectra are still unex- 
plained (see, e.g., Reimers et al [Q in the course of the 
Hamburg ESO survey). 

Investigations on the electronic structure in the pres- 
ence of a magnetic field appear to be quite complicated 
due to the intricate geometry of this quantum prob- 
lem. For the hydrogen atom the impact of the com- 
peting Coulomb and diamagnetic interaction is particu- 
larly evident and pronounced in the intermediate regime 
for which the magnetic and Coulomb forces are compa- 
rable. For different electronic degrees of excitation of 
the atom the intermediate regime is met for different 
absolute values of the field strength. For the ground 
state this regime corresponds to field strengths around 
J — 1 (for the magnetic field strength as well as for 
other physical values we use atomic units and, in partic- 
ular, 7 = B/Bq, Bq corresponds to the magnetic field 
strength Bq = hc/eal = 2.3505-10'5T). Both early |3|jl4| 
and more recent works |Q,0 on the hydrogen atom have 
used different approaches for relatively weak fields (the 
Coulomb force prevails over the magnetic force) and for 
very strong fields (the Coulomb force can be considered 
as weak in comparison with the magnetic forces which is 
the so-called adiabatic regime). A powerful method to 
obtain comprehensive results on low-lying energy levels 
of the hydrogen atom in particular in the intermediate 
regime is provided by mesh methods B . For atoms with 
several electrons there are two decisive factors which en- 
rich the possible changes in the electronic structure with 
varying field strength compared to the one-electron sys- 
tem. First we have a third competing interaction which 
is the electron-electron repulsion and second the differ- 



ent electrons feel very different Coulomb forces, i.e. pos- 
sess different one particle energies, and consequently the 
regime of the intermediate field strengths appears to be 
the sum of the intermediate regimes for the separate elec- 
trons. 

Opposite to the hydrogen atom the wavefunctions of 
the multi-electron atoms change their symmetries with 
increasing field strength. It is well known that the sin- 
glet zero-field ground state of the helium atom (Is^ in 
the Hartree-Fock language) is replaced in the high-field 
regime by the triplet fully spin polarised configuration 
ls2p_i. For atoms with more than two electrons the 
evolution of the ground state within the whole range of 
field strengths < 7 < -f cxo includes multiple intermedi- 
ate configurations besides the zero-field ground state and 
the ground state corresponding to the high field limit. In 
view of the above there is a need for further quantum 
mechanical investigations and data on atoms with more 
than two electrons in order to understand their electronic 
structure in strong magnetic fields. Our calculations al- 
lowed us to obtain the first conclusive results on the series 
of ground state configurations for the Li ||l^ and C ||l7| ] 
atoms. These results are substantially different from pre- 
viously published ones |1^. The ground state electronic 
configurations of the beryllium atom for < "f < +00 
were not investigated so far. A previous work on the 
beryllium atom ||T^ focused on problems associated with 
the symmetries of the Hartree-Fock wavefunction of the 
low-field ground state ls^2s^ of this atom. For strong 
fields the ls^2s^ state represents a highly excited state 
and the electronic ground state configuration of Be is, so 
far, not investigated. 

In the current paper we present results of our fully nu- 
merical 2D Hartree-Fock mesh calculations of the beryl- 
lium atom and Be+ ion in magnetic fields and obtain for 
the first time conclusive results on the structure and en- 
ergy of the ground state configurations of these systems 
for arbitrary field strengths. 

II. METHOD 

The computational method applied in the current 
work coincides with the method described in our works 
ri-pl and applied afterwards in [|l6|,|7|j23|J2l . We 



solve the electronic Schrodinger equation for the beryl- 
lium atom in a magnetic field under the assumption of 
an infinitely heavy nucleus in the (unrestricted) Hartree- 
Fock (HF) approximation. The solution is established in 
the cylindrical coordinate system (p, (j), z) with the z-axis 
oriented along the magnetic field. We prescribe to each 
electron a definite value of the magnetic quantum num- 
ber 771^. Each one-electron wave function ^^ depends on 
the variables (^ and (p, z) 



*^(p,</),z) = (2^)-i/^e-™''*VM(z,p) 



(1) 



where /i indicates the numbering of the electrons. The re- 
sulting partial differential equations for tp^^z^p) and the 
formulae for the Coulomb and exchange potentials have 
been presented in ref. ua. These equations as well as 
the Poisson equations for inter-electronic Coulomb and 
exchange potentials are solved by means of the fully nu- 
merical mesh method described in refs. 0,12^. The finite- 
difference solution of the Poisson equations on sets of 
nodes coinciding with those of the Hartree-Fock equa- 
tions turns out to be possible due to a special form of 
uniform meshes used in the present calculations and in 
refs. |16| , p7| , p^ . Details and discussion on these meshes 
are presented in ref. pq |. 

Our mesh approach is flexible enough to yield precise 
results for arbitrary field strengths. Some minor decrease 
of the precision appears for electronic configurations with 
big differences in the spatial distribution of the electronic 
density for different electrons. This results in big differ- 
ences with respect to the spatial extension of the density 
distribution for different electrons. This situation is more 
typical for the electronic configurations which do not rep- 
resent the ground state at the corresponding fields (e.g. 
ls^2s^ at very strong fields or ls2p_i3(i-24/-3 in the 
weak field regime). The precision of our results depends, 
of course, on the number of mesh nodes and can be al- 
ways improved in calculations with denser meshes. Most 
of the present calculations are carried out on sequences of 
meshes with the maximal number of nodes being 80 x 80. 

Along with the numerical solution of the Schrodinger 
equation the key element for solving the problem of the 
ground state electronic configurations is a proper choice 
of the configurations, which could potentially be the 
ground state ones. An example of solving this problem 
is presented in ||l^]. In that work we have developed a 



strategy which enables one to reduce the set of possi- 
ble ground state configurations which are then subject 
to a following numerical investigation. This removes the 
risk of missing some ground state configurations due to 
the limited possibilities of performing numerical investi- 
gations. With increasing number of electrons the number 
of configurations which cannot a priori be excluded from 
becoming the ground state increase rapidly. A compre- 
hensive numerical investigation of all these configurations 
is, in general, not feasible. The above-mentioned strat- 
egy to exclude certain configurations is therefore highly 
desirable. It is based on a combination of qualitative 
theoretical arguments and numerical calculations of the 
energies of electronic configurations. As a first step the 
set of electronic configurations has to be separated into 
several groups according to their spin projections Sz- 
The following considerations have to be carried out in 
each subset separately, and are certainly more transpar- 
ent by starting with the limit of infinite strong fields and 
analysing the electronic configurations with decreasing 
field strengths. The qualitative theoretical considerations 
mentioned above are based on the geometry of the spa- 
tial part of the wavefunction and enable one to determine 
the ground state for the high-field limit as well as several 
candidates for the ground state configuration with de- 
creasing field strength. The numerical calculations then 
enable us to decide which of these candidates becomes 
the actual first intermediate ground state and yields the 
transition field strength. The knowledge of the first inter- 
mediate ground state allows us to repeat the qualitative 
considerations for the second intermediate ground state 
to obtain a list of candidates which is then investigated 
by means of numerical calculations. Repeating this pro- 
cedure one can determine the full sequence of the ground 
state configurations for each subset ^2 and finally the 
sequence of ground state configurations for the physical 
system. 

III. GROUND STATE ELECTRONIC 
CONFIGURATIONS FOR 7 = AND 7 ^ 00 

In this section we provide some helpful qualitative con- 
siderations on the problem of the atomic multi-electron 
ground states particularly in the limit of strong magnetic 



fields. 

For the case 7 = the ground state configuration of the 
beryllium atom can be characterised as ls^2s^. This no- 
tation has a literal meaning when considering the atom in 
the framework of the restricted Hartree-Fock approach. 
The latter is an approximation of limited quality in de- 
scribing the beryllium atom as it was shown in many fully 
correlated calculations both for the field-free Be atom 
p6| , p7t and for its polarizabilities in electric fields pq,p9i . 
It was pointed out in these works that the Be atom is a 
strongly correlated system and that the HF ground state 
wavefunction (i. e. the spherically symmetric ls^2s^) is 
not a very accurate zeroth-order wavefunction, especially 
for calculations of electric polarizabilities. This is due to 
a significant contribution of the ls^2p^ configuration to 
the ground state wave function. The latter configuration 
is evidently a non-spherical one. This fact is in agree- 
ment with results of ref. p9| where the fully numerical 
2D unrestricted Hartree-Fock approach provides the 2s^ 
shell stretched along the z axis even for 7 = 0. In terms 
of spherical functions it is natural to describe this geom- 
etry of the 2s^ shell as a mixture of 2s and 2pQ functions. 
We remark that the s, p, d. . . orbital notation both for 
7 = and 7 7^ is based on the behaviour of the wave 
functions in the vicinity of the origin and on the topol- 
ogy of the nodal surfaces, but does not imply any detailed 
geometry or certain values of the orbital moment I. 

It is evident that the field-free ground state of the 
beryllium atom remains the ground state only for rela- 
tively weak fields. The set of one-electron wave functions 
constituting the HF ground state for the opposite case 
of extremely strong magnetic fields can be determined as 
follows. The nuclear attraction energies and HF poten- 
tials (which determine the motion along z axis) are small 
for large 7 in comparison to the interaction energies with 
the magnetic field (which determines the motion perpen- 
dicular to the magnetic field and is responsible for the 
Landau zonal structure of the spectrum). Thus, all the 
one-electron wavefunctions must correspond to the low- 
est Landau zones, i.e. the magnetic quantum numbers 
771^ are not postive for all the electrons m^ < 0, and the 
system must be fully spin-polarised, i.e. s^^ = —^- For 
the Coulomb central field the one-electron levels form (as 
B —f 00) quasi ID Coulomb series with the binding en- 



ergy es = 2^ for n^ > and e^ — > 00 for n^ = 0, 
where n^ is the number of nodal surfaces of the wave 
function with respect to the z axis. The binding energy 
of a separate electron has the form 



es = (m + \m\ + 2s^ + 1)7/2 - 



(2) 



where e is the energy of the electron. 

When considering the case 7 — > cxd it is evident, that 
the wave functions with n^ = have to be chosen for the 
ground state configuration. Furthermore starting with 
the energetically lowest one particle level the electrons 
occupy according to the above arguments orbitals with 
increasing absolute value of the magnetic quantum num- 
ber m^. Consequently the ground state of the beryllium 
atom must be given by the fully spin-polarised configu- 
ration ls2p_i3d_24/-3. In our notation of the electronic 
configurations we assume in the following that all paired 
electrons, like for example the Is^ part of a configuration, 
are of course in a spin up and spin down orbital, respec- 
tively, whereas all unpaired electrons possess a negative 
projection of the spin onto the magnetic field direction. 
On a qualitative level the configuration ls2p_i3(i_24/_3 
is not very different from similar electronic configurations 
for other atoms (see ref. |2J|). This is a manifestation of 
the simplification of the picture of atomic properties in 
the limit 7 — > cxd where a linear sequence of electronic 
configurations replaces the periodic table of elements of 
the field-free case. 

The problem of the configurations of the ground state 
for the intermediate field region cannot be solved without 
doing explicit calculations combined with some qualita- 
tive considerations in order to extract the relevant con- 
figurations. 

IV. GROUND STATE ELECTRONIC 

CONFIGURATIONS FOR ARBITRARY FIELD 

STRENGTHS 

In order to determine the ground state electronic con- 
figurations of the beryllium atom wc employ here the 
strategy introduced in ref. |jl^ where the carbon atom 
has been investigated. First of all, we divide the possible 
ground state configurations into three groups according 



to their total spin projection Sz '■ the Sz = group (low- 
field ground state configurations) , the intermediate group 
5"^ = — 1 and the Sz = —2 group (the high- field ground 
state configurations). This grouping is required for the 
following qualitative considerations which are based on 
the geometry of the spatial parts of the one-electron wave 
functions. In the course of this discussion it is expedi- 
ent to treat local ground states for each Sz subset (i.e. 
the lowest states with a certain Sz value) along with the 
global ground state of the atom as a whole. For each value 
of the magnetic field strength one of these local ground 
states is the global ground state of the atom. 

According to the general arguments presented in the 
previous section we know that the ground state configu- 
ration of the beryllium atom in the high field limit must 
be the fully spin-polarised state ls2p_i3(i_24/_3. The 
question of the ground state configurations at interme- 
diate fields cannot be solved without performing explicit 
electronic structure calculations. On the other hand, the 
a priori set of possible intermediate ground state con- 
figurations increases enormously with increasing number 
of electrons and is rather large already for the beryl- 
lium atom. Some qualitative considerations are there- 
fore needed in order to exclude certain configurations as 
possible ground state configurations thereby reducing the 
number of candidates for which explicit calculations have 
to be performed. As mentioned in the previous section 
the optimal strategy hereby consists of the repeated pro- 
cedure of determining neighbouring ground state config- 
urations with increasing (or decreasing) magnetic field 
strength using both qualitative arguments as well as the 
results of the calculations for concrete configurations. 

The total energies for the considered states and par- 
ticularly of those states which become the global ground 
state of the atom for some regime of the field strength 
are illustrated in figure 1 . In the following paragraphs we 
describe our sequence of selection procedure and calcu- 
lations for the candidates of the electronic ground state 
configurations. 

Due to the simplicity of the ground state electronic 
configurations of atoms in the limit 7 —> 00 it is natural 
to start the consideration for 7 y^ with the high-field 
ground state and then consider other possible candidates 
in question for the electronic ground state for S'^ = — 2 



(see figure 1) with decreasing field strength. The con- 
sideration of the high-field (i.e. the fully spin-polarised) 
regime was carried out in ref. [p4| ^ind for this case (i.e. 
Sz = —2 for beryllium) we repeat this consideration in 
more detail. In particular, we have found in ref. |2J] that 
the beryllium atom, opposite to the carbon and heavier 
elements has only one fully spin-polarised ground state 
configuration. 

All the one electron wave functions of the high-field 
ground state ls2p_i3d_24/_3 possess no nodal surfaces 
crossing the z-axis and occupy the energetically lowest 
orbitals with magnetic quantum numbers ranging from 
m ^ down to to = —3. The 4/_3 orbital possesses 
the smallest binding energy of all orbitals constituting 
the high-field ground state. Its binding energy decreases 
rapidly with decreasing field strength. Thus, we can 
expect that the first crossover of ground state config- 
urations happens due to a change of the 4/_3 orbital 
into one possessing a higher binding energy at the cor- 
responding lowered field strength. One may think that 
the first transition while decreasing the magnetic field 
strength will involve a transition from an orbital possess- 
ing n^ = to one for n^ — 1. The energetically lowest 
available one particle state with n^ = 1 is the 2po or- 
bital. Another possible orbital into which the 4/_3 wave 
function could evolve is the 2s state. For the hydro- 
gen atom or hydrogen-like ions in a magnetic field the 
2po is stronger bound than the 2s orbital. On the other 
hand, owing to the electron screening in multi-electron 
atoms in field-free space the 2s orbital tends to be more 
tightly bound than the 2po orbital. Thus, two states i.e. 
ls2po2p-i3(i_2 and Is2s2p_i3(i_2 are candidates for be- 
coming the ground state in the Sz — —2 set when we 
lower the field strength coming from the high field sit- 
uation. The numerical calculations show that the first 
crossover of the Sz ~ —2 subset takes place between the 
ls2p_i3(i_24/-3 and ls2po2p-i3(i_2 configurations (fig- 
ure 1). On the other hand, the calculations show that 
even earlier (i.e. at higher magnetic field strengths) the 
global ground state acquires the total spin S'^ = — 1 due 
to a crossover of the energy curve of the ls2p_i3(i_24/-3 
configuration with that of the configuration ls^2p_i3(i_2 
(which is the local ground state for the Sz — ~^ subset 
in the high- field limit). For the fields below this point 



7 = 4.567 the ground state electronic configurations of 
the beryllium atom belong to the subset Sz — —1- This 
means that the beryllium atom has only one fully spin po- 
larised ground state configuration (as mentioned above). 

The electronic configurations 

ls^2p_i3(i_2 and ls2p_i3(i_24/-3 differ by the replace- 
ment of the spin down 4/_3 orbital through the spin up Is 
orbital and according to the argumentation presented in 
the previous section the ls^2p_i3ci_2 represents the local 
ground state configuration for the subset Sz = — 1 in the 
limit 7 — > cxD. Analogous arguments to that presented in 
the previous paragraph provide the conclusion, that in 
the process of decreasing field strength the ls^2p_i3(i_2 
ground state electronic configuration can be replaced ei- 
ther by the ls^2s2p_i or by the ls^2po2p-i configura- 
tion. The numerical calculations show, that the curve 
^Is22s2j9_i(7) intersects the curve -Bis22p_i3d_2(7) at a 
higher magnetic field (7 = 0.957) than £'is2 2po2p_i(7) 
crosses -Eis22p_i3d_9(7)- The difference with respect to 
the order of the local ground state configurations in the 
subsets Sz = —2 and 5*2 = —1 stems from the differ- 
ence in the magnetic field strengths characteristic for the 
crossovers in these subsets. At moderate field strengths 
{Sz = — 1) the influence of the Coulomb fields of the 
nucleus and electrons prevails over the influence of the 
magnetic field and make the energy of the 2s orbital 
lower than that of the 2j>o orbital. On the other hand, 
at stronger fields characteristic for the subset Sz = —2 
the energies of these orbitals are governed mostly by the 
magnetic field and, in result, the energy of the 2po orbital 
becomes lower than the energy of the 2s orbital. 

From our simple qualitative considerations we can 
conclude, that the configuration ls^2s2p_i is the local 
ground state configuration of the subset 5*2 = —1 for 
the weak field case, i.e. for 7 — > 0. Indeed, when we 
construct such a configuration, the first three electrons 
go to orbitals Is and 2s forming the ls^2s configuration 
with Sz = —1/2. The fourth electron must then have the 
same spin as the 2s orbital electron to obtain the total 
spin value Sz — —1- Thus, the lowest orbital which it can 
occupy is the 2p_i. Therefore, there are two local ground 
state configurations in the subset 5*2 = —1 and they both 
represent the global ground state for some ranges of the 
magnetic field strengths. 



The necessary considerations for the subset Sz = 
are quite simple. At 7 = and, evidently, for very weak 
fields the ground state of the beryllium atom has the con- 
figuration ls^2s^. We can expect, that when increasing 
the magnetic field strength, the next lowest state with 
Sz — will be the ls^2s2p_i configuration with oppo- 
site directions of the spins of the 2s and 2p_i electrons. 
But both contributions, the Zeeman spin term and the 
electronic exchange make the energy of this state higher 
than the energy of the state ls^2s2p_i with the parallel 
orientation of the spins of the 2s and 2p_i electrons (i.e. 
Sz — —1) considered above. The calculated energies for 
these states are presented in figure 1. Thus, the beryl- 
lium atom has one ground state electronic configuration 
ls^2s^ with the total spin ^-projection Sz = 0. This 
state is the global ground state for the magnetic field 
strengths between 7 = and 7 = 0.0412. Above this 
point the ground state configuration is ls^2s2p_i with 

Sz = -1. 

Summarising the results on the ground state config- 
urations of the beryllium atom we can state that de- 
pending on the magnetic field strength this atom has 
four different electronic ground state configurations. For 
< 7 < 0.0412 the ground state configuration coin- 
cides with the field-free ground state configuration ls^2s^ 
which has zero values for the total magnetic quantum 
number M and spin projection S^. The following are two 
ground state configurations with Sz = —1: ls^2s2p_i 



(M 



-1) for 0.0412 < 7 < 0.957 and ls'^2p^i3d^2 



(M = -3) for 0.957 < 7 < 4.567. For 7 > 4.567 
the ground state configuration is ls2p_i3(i_24/-3 with 
Sz = —2 and M = —6. The complete results of the 
investigations of the sequence of the ground state config- 
urations of the Be atom are presented in table || which 
contains the critical values of 7 at which the crossovers 
of different ground state configurations take place. 

The next aim of this section is the corresponding inves- 
tigation of the ground state configurations of the ion Be+ . 
The field-free ground state of this ion corresponds to the 
ls^2s configuration {Sz = —1/2 and M = 0). In the op- 
posite case 7 — > c» the ground state is obviously given by 
the Is2p_i3(i_2 configuration {Sz — —3/2 and M = —3). 
Thus, we need to investigate only two different subsets 
of electronic ground state configurations: Sz = —1/2 and 



^2 = —3/2. The energy curves which are necessary for 
this investigation are presented in figure 2. The subset 
Sz = —1/2 contains only two possible ground state con- 
figurations ls^2s and ls^2p_i. The latter is the local 
ground state configuration for this subset in the limit 
7 —^ 00. The curves -Eis22s(7) f^nd Eig22p_-^{^) intersect 
at 7 = 0.3185 and above this point £'is22p_i < £'is22s- In 
the subset Sz — —3/2 we have to consider the configura- 
tions ls2po2p-i and Is2s2p_i along with the high-field 
ground state configuration Is2p_i3(i_2. But the numeri- 
cal calculations show that the energies of both ls2po2p-i 
and Is2s2p_i lie above the energy of the Is2p_i3(i_2 con- 
figuration at the intersection point (7 — 4.501) between 
Eis2p-i3d-2il) and i?is22p_i(7)- Thus, the ion Be+ has 
three different electronic ground state configurations in 
external magnetic fields: for < 7 < 0.3185 it is ls^2s 
{Sz = -1/2 and A/ = 0), then for 0.38483 < 7 < 4.501 
it is ls22p_i {Sz = -1/2 and M = -1) and for all 
the values 7 > 4.501 the ground state configuration is 
Is2p_i3d_2 {Sz = -3/2 and M = -3). These results 
are summarised in table y. The set of the electronic 
ground state configurations for the Be+ ion appears to be 
qualitatively the same as for the lithium atom ||l^. The 
field strengths for the corresponding transition points are 
roughly two times higher for the Be^ ion than for the Li 
atom. 

V. SELECTED QUANTITATIVE ASPECTS 



In tables III and IV we present the total energies of the 
four ground state electronic configurations of the beryl- 
lium atom and the three ground state electronic config- 
urations of the ion Bc+, respectively. These data cover 
a very broad range of the field strengths from 7 = and 
very weak magnetic fields starting with 7 = 0.001 up to 
extremely strong fields 7 = 10000. The latter value of 
the field strength can be considered as a rough limit of 
applicability of the non-rclativistic quantum equations to 
the problem (see below) . The corresponding data on the 
Be+ ion can be found in tables || and IV . 

So far there exist three works which should be men- 
tioned in the context of the problem of the beryllium 
atom in strong magnetic fields. Rcf. 11^] deals with the 
ls^2s^ state of this atom in fields < 7 < 1000 and 



ref. pj] investigates the ground state energies of atoms 
with nuclear charge number Z < 10 in the high-field, i.e. 
fully spin polarised regime. Both these works contain 
calculations carried out by the method used in the cur- 
rent work and do not represent a basis for comparison. 
The comparison of our results with an adiabatic Hartree- 
Fock calculation of atoms with Z < 10 ^^ is presented 
in pM and we can briefly summarise this comparison for 
two values of the magnetic field strengths: for B12 — 0.1 
(i.e. B ^ O.lx IQi^G) our result is £; = -0.89833keV 
whereas ref. @ yields E = -0.846keV; for Bu = 5 (i.e. 
B = 5x lO^^Q) ouj. j-esuh is £: = -3.61033keV, whereas 
ref. |32l yields E = — 3.5840keV. This comparison allows 
us to draw the conclusion of a relatively low precision 
of the adiabatic approximation for multi-electron atoms 
even for relatively high magnetic fields. 

In figure 3 we present the ionization energy i?ion of the 
beryllium atom depending on the magnetic field strength. 
This continuous dependence is divided into six parts cor- 
responding to different pairs of the ground state config- 
urations of the Be atom and Be"*" ion involved into the 
ionization energy. The five vertical dotted lines in figure 3 
mark the boundaries of these sections. The alteration of 
the sections of growing and decreasing ionization energy 
originates from diiferent dependencies of the total ener- 
gies of the Be and Be+ on the magnetic field strength for 
different pairs of the ground state configurations of these 
two systems. One can see the sharp decrease of the ion- 
ization energy between the crossovers (4) and (5). This 
behaviour is due to the fact that i^ion is determined in 
this section by the rapidly decreasing total energy of the 
state Is2p_i3(i_2 of the Be"*" ion (figure 2) and by the 
energy of the Be atom in the state ls^2p_i3fi_2 which 
is very weakly dependent on the field strength (figure 
1). Another remarkable feature of the curve -Eion(7) is 
its behaviour in the range of field strengths between the 
transitions (2) and (3). The ionization energy in this re- 
gion contains a very shallow maximum and in the whole 
section it is almost independent on the magnetic field. 
Thus, the ionization energy is stationary in this regime 
of field strengths 7 = 0.3 — 0.5 a.u. typical for many 
magnetic white dwarfs |^. 

The above-discussed properties are based on the be- 
haviour of the total energy of the Be atom and Be+ ion. 



On the other hand, the behaviour of the wavefunctions 
and many intrinsic characteristics of atoms in external 
magnetic fields are associated not with the total energy, 
but with the binding energies of separate electrons (g) 
and the total binding energy of the system 



N 

E 



(mp 



TO. 



2s, 



1)7/2 -£; 



(3) 



where N is the number of electrons. The binding en- 
ergies of the ground state electronic configurations of 
Be and Be+ depending on the magnetic field strength 
are presented in figures 4 and 5. These dependencies 
at very strong magnetic fields may illustrate our consid- 
erations of the previous sections. One can see in fig- 
ure 4 that the high-field ground state ls2p_i3d_24/_3 
is not the most tightly bound state of the beryllium 
atom. For all the values of the magnetic fields consid- 
ered in this paper its binding energy is lower than that 
of states ls^2s2p_i and ls^2p_i3d_2 and for 7 < 100 it 
is lower than E^ig22s2. The latter circumstance can be 
easily explained by the fact that the ls^2s^ configura- 
tion contains two tightly bound orbitals Is whereas the 
ls2p_i3d_24/-3 possess only one such orbital. However, 
with increasing magnetic field strengths the contribution 
of the group 2p_i3(i_24/-3 to the binding energy turns 
out to be larger than that of the ls2s^ group. Analo- 
gously we can expect SBis2p_i3rf_24/_3 > ■E^Bis22s2p_i at 
some very large fields 7 > 10000. On the other hand, 
it is evident that the state ls2p_i3(i_24/-3 must be 
less bound than ls^2p^i3d^2 because both these con- 
figurations are constructed of orbitals with binding en- 
ergies, logarithmically increasing as 7 — > 00, but the 
ls^2p_i3(i_2 contains an additional Is orbital, which is 
more tightly bound than 4/_3 at arbitrary field strengths. 
The plot for the Be+ ion (figure 5) illustrates the same 
features and one can see in this figure nearly parallel 
curves Sbis22p_i(7) and -EBis2p_i3d_2(7) in the strong 
field regime. 

Figures 6 and 7 allow us to add some details to the 
considerations of the previous section. These figures 
present spatial distributions of the total electronic densi- 
ties for the ground state configurations of the beryllium 
atom and its positive ion, respectively. These pictures 
allow us to gain insights into the geometry of the dis- 
tribution of the electronic density in space and in par- 



ticular its dependence on the magnetic quantum num- 
ber and the total spin. The first pictures in these fig- 
ures present the distribution of the electronic density 
for the ground state for 7 = 0. The following pictures 
show the distributions of the electronic densities for val- 
ues of the field strength which mark the boundaries of 
the regimes of field strengths belonging to the different 
ground state configurations. For the high-field ground 
states we present the distribution of the electronic density 
at the crossover field strength and at 7 = 500. For each 
configuration the effect of the increasing field strength 
consists in compressing the electronic distribution to- 
wards the z axis. However the crossovers of ground state 
configurations involve the opposite effect due to the fact 
that these crossovers are associated with an increase of 
the total magnetic quantum number M . 

In the first rows of figures 6 and 7 one can see a dense 
core of Is^ electrons inside the bold solid line contour 
and a diffuse distribution of 2 s electrons outside this 
core. The prolate shape of the bold solid line contour 
in the first plot of the figure 6 (ls^2s^, 7 = 0) reflects 
the non-spherical distribution of the 2s electrons in our 
calculations or the admixture of the ls^2pQ configuration 
to the ls^2s^ one from the point of view of the multi- 
configurational approach [p6|-p9[. 

Some additional issues concerning the results presented 
above have to be discussed. First, our HF results do 
not include the effects of correlation. To take into ac- 
count the latter would require a multi-configurational ap- 
proach which goes beyond the scope of the present paper. 
We, however, do not expect that the correlation energy 
changes our main conclusions like, for example, the fact of 
the crossovers with respect to the different ground states 
configurations. With increasing field strength the effec- 
tive one particle picture should be an increasingly better 
description of the wave function and the percentage of 
the correlation energy should therefore decrease (see ref. 
p3[ for an investigation on this subject). Two other im- 
portant issues are relativistic effects and effects due to 
the finite nuclear mass. Both these points are basically 
important for very high magnetic field strengths and they 
have been discussed in ref. |g^. For the systems Be and 
Be+ and for most of the field strengths considered here 
these effects result in minor corrections to the energy. 



VI. SUMMARY AND CONCLUSIONS 

We have applied our 2D mesh Hartree-Fock method 
to the magnetised neutral beryllium atom and beryllium 
positive ion. The method is fiexible enough to yield pre- 
cise results for arbitrary field strengths and our calcula- 
tions for the ground and several excited states are per- 
formed for magnetic field strengths ranging from zero up 
to 2.3505 • lO^T (7 = 10000). Our considerations focused 
on the ground states and their crossovers with increasing 
field strength. The ground state of the beryllium atom 
undergoes three transitions involving four different elec- 
tronic configurations. For weak fields up to 7 = 0.0412 
the ground state arises from the field-free ground state 
configuration ls^2s^ with the total spin z-projection 
Sz = 0. With increasing strength of the field two differ- 
ent electronic configurations with Sz — ~^ consequently 
become the ground state: ls^2s2p_i and ls^2p_i3(i_2- 
At 7 = 4.567 the last crossover of the ground state con- 
figurations takes place and for 7 > 4.567 the ground state 
wavefunction is represented by the high-field-limit fully 
spin polarised configuration ls2p_i3(i_24/_3, Sz ~ —2. 

For the ion Be^ we obtain three different ground state 
configurations possessing two values of the spin projec- 
tion. For fields below 7 = 0.3185 the ground state 
electronic configuration has the spin projection Sz — 
— 1/2, magnetic quantum number M = and qualita- 
tively coincides with the zero-field ground state config- 
uration ls^2s. Between 7 = 0.3185 and 7 = 4.501 the 
ground state is represented by another configuration with 
Sz = -1/2, i.e. ls'^2p^i {M = -1). Above the point 
7 — 4.501 the fully spin polarised high-field-limit config- 
uration Is2p_i3d_2 {Sz = —3/2) is the actual ground 
state of the Be+ ion. Thus, the sequence of electronic 
ground state configurations for the Be"*" ion is similar to 
the sequence for the Li atom [|6|. We present detailed 
tables of energies of the ground state configurations for 
Be and Be"*". 

For Be and Be+ we have presented also the binding 
energies of the ground state configurations dependent on 
the magnetic field strength and maps of electronic densi- 
ties for these configurations. For the Be atom we present 
its ionization energy dependent on the field strength. 

Our investigation represents the first conclusive study 



of the ground state of the berylhum atom and Be+ ion 
for arbitrary field strengths. For the Be atom we have 
obtained a new sequence of electronic configurations with 
increasing field strength. This sequence does not coincide 
with any such sequences obtained previously for other 
atoms and ions and could not be predicted even quali- 
tatively without detailed calculations. Putting together 
what we currently know about ground states of atomic 
systems in strong magnetic fields we can conclude that 
the H, He, Li, Be, C, He+, Li+ and Be"*" ground states 
have been identified. For other atoms and multiple series 
of ions the question about the ground state configurations 
is still open. 
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Figure Captions 

Figure 1. The total energies (in atomic units) of the 
states of the berylhum atom as functions of the magnetic 
field strength considered for the determination of the 
ground state electronic configurations. The field strength 
is given in imits of 7 = (-|;),Bo = fic/eal = 2.3505-10''^T. 

Figure 2. The total energies (in atomic units) of the 
states of the beryllium positive ion as functions of the 
magnetic field strength considered for the determination 
of the ground state electronic configurations. The field 
strength is given in units of 7 = {^),Bq = hc/euQ = 
2.3505-lO^T. 

Figure 3. Be atom ground state ionization energy 
Ei. Transition points are marked by broken vertical 
lines. The sequence of the transitions are (from left to 
right): 1. Be: ls'^2s'^ — > ls'^2s2p^i; 2. Bc+: ls'^2s — > 
ls22p_i; 3. Be: ls'^2s2p^i — > ls'^2p^i3d^2- 4. Bc+: 
ls22p_i — > Is2p_i3d_2; 5. Be: ls'^2p^i3d^2 — > 
ls2p_i3(i_24/-3. Crossovers (4) and (5) take place at 
relatively close values of 7 and are not resolved in the 
figure. 

Figure 4. The binding energies (in atomic units) 
of the ground state electronic configurations of the Be 
atom depending on the magnetic field strength. The field 
strength is given in units of 7 = (-g-), -Bo = Tic/ea^ — 
2.3505-lO^T. 

Figure 5. The binding energies (in atomic units) 
of the ground state electronic configurations of the Be"*" 
ion depending on the magnetic field strength. The field 
strength is given in units of 7 = (-^), -Bo = Tic/eaf^ = 
2.3505-lO^T. 

Figure 6. Contour plots of the total electronic den- 
sities for the ground state of the beryllium atom. For 
neighbouring lines the densities are different by a factor of 
2. The coordinates z, x as well as the corresponding field 
strengths are given in atomic units. Each row presents 
plots for a ground state configuration at its lower (left) 
and upper (right) intersection points. Rows: 1. ls^2s^: 
7 = and 7 = 0.0412; 2. Is'^2s2p^i: 7 = 0.0412 and 
7 = 0.957; 3. ls^2p^iM-2-- 7 = 0.957 and 7 = 4.567; 4. 



ls2p_i3<i_24./-3: 7 == 4.567 and 7 = 500. 

Figure 7. Contour plots of the total electronic den- 
sities for the ground state of the beryllium positive ion. 
For neighbouring lines the densities are different by a 
factor of 2. The coordinates z, a; as well as the corre- 
sponding field strengths arc given in atomic units. Each 
row presents plots for a ground state configuration at its 
lower (left) and upper (right) intersection points. Rows: 
1. ls^2s: 7 == and 7 = 0.3185; 2. ls'^2p^i: 7 = 0.3185 
and 7 = 4.501; 3. ls2p^iM-2- 7 = 4.501 and 7 = 500. 
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TABLE I. The Hartree-Fock ground state configurations of tlie beryllium atom in external magnetic fields. The configurations 
presented in the table are the ground state configurations for 7min < 7 < 7max- Atomic units are used. 



7min 



7niax 



The ground state configuration 



M 



S, 



-E(7min) 



1 





0.0412 


1 2o 2 

Is 2s 


2 


0.0412 


0.957 


ls^2s2p_i 


3 


0.957 


4.567 


ls^2p-iM-2 


4 


4.567 


c» 


ls2p_i3d-24/-3 









-14.57336 


-1 


-1 


-14.57098 


-3 


-1 


-15.13756 


-6 


-2 


-15.91660 



TABLE IL The Hartree-Fock ground state configurations of the Be"*" ion in external magnetic fields. The configurations 
presented in the table are the ground state configurations for 7min < 7 < 7max. Atomic units are used. 



The ground state configuration 



M 



S, 



E{jmin) 



1 





0.3185 


ls^2s 


2 


0.3185 


4.501 


ls^2p 


3 


4.501 


oo 


ls2p- 






-1/2 


-14.27747 


1 


-1/2 


-14.38602 


3 


-3/2 


-15.01775 
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TABLE III. The total energies of 
strength. The figures in parentheses 
Atomic units are used. 



the ground state configurations of the beryllium atom depending on the magnetic field 
are the labels of the ground state configurations provided in the first column of table I. 



7 



E{1) 



E{2) 



m) 



S(4) 



0.000 

0.001 

0.002 

0.005 

0.01 

0.02 

0.03 

0.04 

0.05 

0.07 

0.1 

0.12 

0.15 

0.2 

0.3 

0.3185 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1. 

1.2 

1.5 

2. 

3. 

4. 

4.501 

5. 

7. 

8. 

10. 

12. 

15. 

20. 

30. 

40. 

50. 

100. 

200. 

500. 



-14.57336 
-14.57336 
-14.57335 
-14.57332 
-14.57322 
-14.57279 
-14.57209 
-14.57111 
-14.56986 
-14.56657 
-14.55971 
-14.55395 
-14.54367 
-14.52261 
-14.46861 

-14.40279 
-14.32860 
-14.24832 
-14.16352 
-14.07526 
-13.98431 
-13.89120 
-13.69990 
-13.40329 
-12.88908 
-11.79811 
-10.633617 

-9.40602 

-6.79760 

-5.43095 

-2.5988 

+0.34064 

+4.9055 

+12.8201 

+29.3964 

+46.5935 

+64.186 

+155.286 

+343.899 

+924.20 



-14.51206 
-14.51357 
-14.51507 
-14.51953 
-14.52690 
-14.54138 
-14.55553 
-14.56933 
-14.58281 
-14.60879 
-14.64548 
-14.66851 
-14.70108 
-14.75065 
-14.83520 
-14.84905 
-14.90464 
-14.96264 
-15.01171 
-15.05368 
-15.08989 
-15.12138 
-15.14899 
-15.19498 
-15.24757 
-15.30815 
-15.36376 
-15.34275 

-15.25183 

-14.89530 

-14.64516 

-14.03046 

-13.29115 

-12.00063 

-9.49118 

-3.59324 

+3.04026 

+10.1472 

+49.4177 

+135.659 

+411.830 



-14.19023 

-14.1928 

-14.1952 

-14.2025 

-14.2142 

-14.2361 

-14.2566 

-14.27587 

-14.29437 

-14.32933 

-14.3780 

-14.40838 

-14.45145 

-14.51761 

-14.63369 

-14.73396 

-14.82272 

-14.90262 

-14.97542 

-15.04232 

-15.10422 

-15.16178 

-15.26583 

-15.39926 

-15.57496 

-15.79985 

-15.90161 

-15.91626 

-15.91027 

-15.71644 

-15.53623 

-15.04644 

-14.41743 

-13.27286 

-10.97100 

-5.40704 

+0.95677 

+7.83395 

+46.25962 

+131.4188 

+405.7027 



-9.44321 

-9.4483 

-9.4532 

-9.4675 

-9.4903 

-9.5331 

-9.5735 

-9.6121 

-9.6493 

-9.7198 

-9.82 

-9.8805 

-9.9692 

-10.1081 

-10.36220 

-10.59384 
-10.80901 
-11.01121 
-11.20281 
-11.38545 
-11.56045 
-11.72880 
-12.04863 
-12.49432 
-13.16961 
-14.35016 
-15.38050 

-16.30690 
-17.94005 
-18.67389 
-20.01753 
-21.23057 
-22.86513 
-25.23250 
-29.11102 
-32.28415 
-35.00768 
-45.10519 
-58.08264 
-80.67357 
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1000. 


+1905.14 


2000. 


+3881.5 


5000. 




10000. 





+888.70 
+1860.40 
+4813.56 
+9770.37 



+880.706 
+1850.052 
+4799.35 
+9752.24 



-102.75480 
-129.9790 
-175.2704 
-217.695 
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TABLE IV. 
The figures in 
units are used 



The total energies of the ground state configurations of the Be"*" ion depending on the magnetic field strength, 
parentheses are the labels of the ground state configurations provided in the first column of table II. Atomic 



7 



E(l) 



E{2) 



Ei3) 



0.000 

0.001 

0.002 

0.005 

0.01 

0.02 

0.03 

0.04 

0.0412 

0.05 

0.07 

0.1 

0.12 

0.15 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.957 

1. 

1.2 

1.5 

2. 

3. 

4. 

4.567 

5. 

7. 

8. 

10. 

12. 

15. 

20. 

30. 

40. 

50. 

100. 

200. 



-14.27747 
-14.27797 
-14.27846 
-14.27995 
-14.28241 
-14.28725 
-14.29198 
-14.29659 
-14.29714 
-14.30111 
-14.30981 
-14.32207 
-14.32972 
-14.34047 
-14.35648 
-14.38212 
-14.40046 
-14.41282 
-14.42022 
-14.42350 
-14.42335 
-14.42029 

-14.41478 
-14.39782 
-14.36143 
-14.28225 
-14.08247 
-13.83797 

-13.55019 

-12.85647 

-12.45821 

-11.57652 

-10.59993 

-8.99386 

-6.03364 

+0.59244 

+7.81557 

+15.4261 

+56.5516 

+145.1649 



-14.13093 
-14.13195 
-14.13294 
-14.13593 
-14.14087 
-14.15066 
-14.16030 
-14.16980 

-14.17916 
-14.19746 
-14.22390 
-14.24088 
-14.26542 
-14.30406 
-14.37402 
-14.43599 
-14.49163 
-14.54210 
-14.58821 
-14.63059 
-14.66971 
-14.69069 
-14.70591 
-14.77070 
-14.85169 
-14.95181 
-15.05201 
-15.05004 

-14.96820 

-14.61928 

-14.37080 

-13.75773 

-13.01900 

-11.72840 

-9.217910 

-3.31723 

+3.31895 

+10.42836 

+49.70820 

+135.95916 



-9.41056 
-9.41358 
-9.41657 
-9.42551 
-9.44028 
-9.46939 
-9.49791 
-9.52587 

-9.55332 

-9.60670 

-9.68356 

-9.73294 

-9.80463 

-9.91878 

-10.13144 

-10.32817 

-10.51259 

-10.68705 

-10.85323 

-11.01236 

-11.16540 

-11.31312 

-11.59490 

-11.98978 

-12.59206 

-13.65352 

-14.58615 

-15.07310 

-15.42817 

-16.91814 

-17.58931 

-18.820184 

-19.93310 

-21.43461 

-23.612005 

-27.18373 

-30.10832 

-32.61959 

-41.93414 

-53.90638 
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500. +425.471 +412.14745 -74.73619 

1000. +906.37 +889.0264 -95.07513 

2000. +1883.08 +1860.7100 -120.11947 

5000. +4844.6 +4814.005 -161.7052 

10000. +9809.3 +9770.643 -200.5709 
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